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1 General 


1.1 Background 

Upper Trishuli-1 HEP is constructed to utilize the hydro potential of 
Trishuli River near Dunche and Haku. The construction site is 
situated in a rugged mountain between Tibet in China and the 
Himalaya Mountain, about 70 km away from capital of Nepal, 
Kathmandu, towards the northwest. The project is a run-of-river 
development scheme. The headwork site is at about 280m 
downstream from the confluence of Bhotekoshi and Trishuli rivers. 
The powerhouse site is 500m upstream from Mailun Dobhan where 
Haku VDC is placed. The installed capacity of Upper Trishuli-1 HEP is 
216 MW. It is a year-round power plant capable of generating about 
1,533.4 GWh electricity yearly on average. 

The Adit No.3 is located on the left bank of headrace tunnel, and the 
headrace tunnel chainage at intersection area is H+6562.747m. The 
cross section of headrace tunnel is horseshoe shape and 
4.2mx5.6m size. The elevation of the inlet is EL. 1194.28m, and the 
elevation of outlet is EL.1212.36m. The length of Adit No.3 is 
309.575m and longitudinal slope is 5.86%. 

According to the Contractor’s interpretation and preliminary 
geological survey data, along the tunnel axis, there are interbedded 
mica schist and quartz schist are distributed and they are classified 
as Type IV, Ill,ll as per Q-method from inlet to outlet, the buried depth 
from inlet to outlet is about 50m to 300m, averaged 200m 
approximately. It is speculated that the whole tunnel is mostly dry, 
individual joints surface is wet. No potential continuous sliding 
surface is found on the tunnel slope. 

The longitudinal profile of Adit No.3 is shown in Figure 1.1. The 
preliminary excavation cross section and preliminary support 


pattern are shown in Figure 1.2~1.7. 
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SCALE: 1:500 


LONGITUDINAL PROFILE OF ADIT No.3 


Figure 1.1 Longitudinal Profile Section of Adit No.3 
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Figure 1.2 Typical section of adit No.3 (rock class I) 
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Figure 1.3 Typical section of adit No.3 (rock class II) 
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Figure 1.4 Typical section of adit No.3 (rock class III) 


Page 3 of 57 


DOOSAN 


UPPER TRISHULI-1 HEP 
(216MW) € Dee 


Surrounding Rock Stability Calculation of Adit No.3 
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Figure 1.5 Typical section of adit No.3 (rock class IV) 
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Figure 1.6 Typical section of turning adit (rock class III) 
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Figure 1.7 Typical section of Junction (rock class II) 


1.2 Design Objective 

This calculation only includes rock class I~IV, and the rock class V 
shall be considered separately according to the actual situation. 

This calculation is provided for detailed design stage and aim to 
confirm the objectives which are as following: 

(1) to generally confirm the stability of the tunnel during 
construction; 

(2) to confirm (or adjust) the previously suggested support 


measures in the construction drawings. 


1.3 Purpose and Contents 

(1) Contract documents, Employer’s Requirement; 
(2) Tunnels and shafts in Rock (EM 1110-2-2901); 
(3) Rock Foundations (EM 1110-1-2908); 
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(4) Building Code Requirements for Structural Concrete and 
Commentary (ACI 318) 

(5) Vlachopoulos, N. and Diederichs, M.S. (2009). Improved 
longitudinal displacement profiles for convergence-confinement 
analysis of deep tunnels. Rock Mechanics and Rock Engineering, 
42(2), 131-146; 

(6) The 2008 Kersten Lecture Integration of geotechnical and 
structural design in tunneling; To be presented at the opening 
keynote address, by Dr Evert Hoek, at the University of Minnesota 
56" Annual Geotechnical Engineering Conference to be held in 
Minneapolis on 29 February 2008; 

(7) Handbook Using the Q-system, Rock mass classification and 
support design; 

(8) Paper: Rock-Support Interaction analysis for tunnels in weak rock 
masses; 

(9) https://www.rocscience.com/help/rs2/turorials/rs2 tunnel lining 
design.htm; 

(10) Drawing: Excavation and initial support drawing of adit No.3 
(Drawing No. UT1-C-150-CVL-DG-43004). 


2 Theoretical Background of Calculations 


2.1 Software 

Software Phase 2D will be applied in the calculation. Phase2 is a 
powerful 2D finite element program for soil and rock applications. 
Phase2 can be used for a wide range of engineering projects 
including excavation design, slope stability, groundwater seepage, 
probabilistic analysis, consolidation, and dynamic analysis 


capabilities. 
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Software UNWEDGE is an interactive software for the stability 
analysis of three-dimensional wedge formed by structural 
discontinuity and underground excavation, which is used to analyze 
the underground excavation problem with discontinuous structural 
plane in rock mass. UNWEDGE calculates the safety factor of 
potentially unstable wedge, and can analyze the influence of 


support system on wedge stability. 


2.2 Modelling Conditions, Material Law, Constitutive Model 

The modelling conditions are plane strain. The material law is 
elasto-plastic combined with strain-softening (for rock). The rock 
mass is modelled as a continuum, homogeneous and isotropic. The 
failure criteria is the Mohr-Coulomb failure criteria. What needs to be 
emphasized is that such continuum model with homogeneous and 
isotropic material is a model, only, and does not necessarily reflect 
reality. In reality, the displacements occur in blocky rock masses 
always along discontinuities (joints, foliation etc.), which leads to 
loosening of the rock mass. The determined plastic zone in the 
continuums model represents an area of the rock mass where rock 
blocks were displaced to each other (= loosening of rock mass). In 
reality, it does not represent an area where the material was subject 


to plastic flow, as in soil. 


2.3 Calculation Methodology 

Convergence-confinement analysis for tunnelling is a standard 
approach for preliminary analysis of anticipated wall deformation 
and support design in squeezing ground. 


Convergence-confinement analysis is a widely used tool for 
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preliminary assessment of squeezing potential and support 
requirements for circular tunnels in variety of geological conditions 
and stress states. An analytical plasticity solution such as that 
developed by Carranza-Torres and Fairhurst (2000) is applied to a 
circular opening in an isotropic stress field. An internal pressure, 
initially equal to the in-situ stress is applied on the inside of the 
excavation boundary. The pressure is incrementally relaxed until the 
excavation boundary condition is that of zero normal stress. The 
extent of plastic yielding and thereby, the boundary deformation is 
calculated at each stage of the process. The result is a continuous 
representation of the deformation-internal pressure relationship for 
the tunnel given a particular material strength, deformability, 
dilation and stress state. 

The internal press is, of course not a representation of reality but 
rather a surrogate for the effect of the gradual reduction of the 
radial resistance provided by the initially present tunnel core 
(material inside the tunnel boundary) transitioning to an exposed 
boundary with a progressively distant tunnel face and ultimately a 
long open tunnel with plane strain conditions. The internal pressure 
that is coupled with a given boundary displacement is a measure of 
the amount of support resistance required to prevent further 
displacement at that point in progressive tunnelling model. 

Elastic deformation of the rock mass starts about two diameters 
ahead of the advancing face and reaches its maximum value at 
about two diameters behind the face. At the face position about one 
third of the total radial closure of the tunnel has already occurred 


and the tunnel face deforms inwards as illustrated in Fig. 2.1. 
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Whether or not these deformations induce stability problems in the 
tunnel depends upon the ratio of rock mass strength to the in-site 
stress level. The tunnel displacement profiles for the roof and the 


invert of a tunnel are shown schematically in the Fig. 2.2. 
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Figure 2.1 The shape of the deformed tunnel profile and 


displacement vectors 
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Longitudinal Displacement Profile __ 


^ Radial displacement 
reaches its final value 
at about two tunnel 
diameters behind the 
face 


Inward deformation 
of tunnel face 


— — Radial displacement reaches 
about one third of its final 
value at the tunnel face 


d 


Direction of 
face advance 


—— Displacements start about two 
tunnel diameters ahead of the face 


Figure 2.2 pattern of elastic deformation in the rock mass 
surrounding an advancing tunnel 

The equations can be found in the Kersten Lecture, appendix 1(Hoek 
et.al.,2008) 

According to the formula listed in geological rock classification and 
documentation “Integration of geotechnical and structural design in 
tunneling (Hoek, E., Carranza-Torres, C., Diederichs, M.S. and 
Corkum, 2008)”. The calculation formula of supporting time is as 


follow: 


Ura 3 
P-R /R 
p f 
tot 
Ri 
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Sg" whenX <0 
p _ H max 
= -3d, 
max a 2P 
1-(1- —J)e** whenX > 0 
where, 
R 


r —The maximum plastic zone depth, m; 
"^ The radial displacement of tunnel face, m; 


R, —Tunnel radius, m; 


X —Distance from the initial support to the tunnel face, m; 


“nex —The maximum radial displacement, m. 


With these analytical equations it is possible to estimate the amount 
of displacements at the specified location of support installation, if 
the plastic radius and displacement far from the tunnel face are 
known. Therefore, to get the displacements and the radius of the 
plastic zone far away from the tunnel face (where the internal 
pressure is zero). It is necessary to calculate numerically with Phase 
2 the tunnel without support, first. 

With the determined displacement profile along the tunnel and a 
given location of support installation, the modeler selects the 
corresponding deformation. 

In the next step the tunnel advance is modelled in Phase 2 by 
reducing the internal pressure in the future tunnel until the 
corresponding deformation is reached, when support is installed. 
The modeler installs now the support in the model. The internal 


pressure is then reduced to zero by factor the field stress zero. 
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Further numerical calculation leads to loading of the support. 


2.4 Maximum Unsupported Span 

The maximum unsupported span is calculated by the following 
formula as per EM 2901: X=2 ESR Q^0.4 

Where, X is maximum unsupported span, ESR is excavation support 
ratio, taken as 1.3. Q values and calculated and applied maximum 
unsupported spans are shown Table 2.1 : 


Table 2.1 Maximum unsupported span of tunnel excavation 


Calculated maximum Maximum unsupported 
Rock class | Q Value 
unsupported span(m) | span for application, X(m) 


| 40 11.4 11 
ll 25 9.4 9 
Ill 7 5.7 5 
IV 2.5 3,8 3 


2.5 Ground Characteristic Curve 
Define the relationship between stress release ratio A; and the 


tunnel deformation ui for an advancing tunnel in stress field. A plot 
of ui; versus A; , the curve is based on the assumption that the rock 
at the tunnel face provides an initial support pressure equal to the 
in-situ stress. As the tunnel face advances and the face moves away 
from the section under consideration, the support pressure 
gradually decreases until it reaches zero at some distance behind 


the face. The characteristic curve of class I~IV are shown as below. 
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Characteristic curve of Class I 
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Figure 2.4 Characteristic curve of class Il 
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Characteristic curve of Class III 


Total Displacement ui [m] 
o 
HM 
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Figure 2.5 Characteristic curve of class III 


Characteristic curve of Class IV 


Total Displacement ui [m] 
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Figure 2.6 Characteristic curve of class IV 
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2.6 Longitudinal displacement profile (LDP) for tunnel 

The Longitudinal Displacement Profile (LDP) is required in order to 
establish the relative position of the tunnel face and the section 
under consideration. 

According to the LDP equations, presented above, the LDP curve of 


every rock mass class are shown in the table below. 


LDP Curve of Class | 
0.60 


0.50 


0.40 


0.30 


0.20 


The displacement of vertex (mm) 


0.00 
-15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 


The distance from tunnel advancing face (m) 


Figure 2.7 Longitudinal Displacement Profile curve of class I 
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LDP Curve of Class II 
1.60 
1.40 


1.20 


The displacement of vertex (mm) 


-15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 


The distance from tunnel advancing face (m) 


Figure 2.8 Longitudinal Displacement Profile curve of class Il 


LDP Curve of Class III 
14.00 


12.00 


10.00 


8.00 


6.00 


4.00 


The displacement of vertex (mm) 


0.00 
-15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 


The distance from tunnel advancing face (m) 


Figure 2.9 Longitudinal Displacement Profile curve of class IlI 
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LDP Curve of Class IV 
60.00 


50.00 


40.00 


30.00 


20.00 


The displacement of vertex (mm) 


0.00 
-15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 


The distance from tunnel advancing face (m) 


Figure 2.10 Longitudinal Displacement Profile curve of class IV 


According to Characteristic Curve and LDP Curve of class l-IV, it is 
determined to adopt the relaxation of stress at the supporting time shown as 


below. 


Table 2.2 Relaxation of rock mass 


Plastic Vertex u Stress 
Supportin x 
Rock zone displacemen met Relaxation 
g time X 
Class radius Rp |t Closure | Percentage 
(m) 
(m) Ug (mm). | ralo (96) 
| 11 3.62 0.5 0.99 95 
II 9 4.44 1.5 0.97 95 
Hl 5 4.45 11.8 0.86 90 
IV 3 8.84 62.5 0.53 85 


Relaxation refer to that (100-Relaxation)% of surrounding rock stress is 
undertaken by rock and support, Relaxation 96 release during excavation. The 


stress relaxation percentage has already considered the margin of safety. 
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3 Properties of Materials 


3.1 Rock Mass Properties 
The mechanical parameters of surrounding rock are shown in Table 
3.3. 


Table 3.1 Mechanical parameter of surrounding rock 


Uniaxial Rock mass/rock a 
: , Uniaxial 
compressi | Elastic l mass . 
Class of Poisson' . compressi | Tensile 
. . ve modulu . shear (shearing) 
surroundin| Unit s ratio ve Strength 
. strength S strength 
g weight of rock strength of! of rock 
(UCS) | Of rock 
rock mass rock mass,| mass 
of intact | mass Peak Residual 
(UCS) 
rock 
Y 
Oc En D peak Cpeak Qa Cres c Oo 
kN/ MPa MPa A m 
(MPa) (GPa) E (°) ( ) |e) ( ) (MPa) (MPa) 
m?) 
| 27.0 60 30.0 0.20 | 50 | 5.0 | 40 |3.33 27.5 2.75 
ll 27.0 60 20.5 0.25 | 48 | 2.5 | 38 |1.67 13.0 1.3 
llla 26.5 60 10.0 0.25 | 45| 1 | 36 |0.66 4.8 0.48 
Ill b 26.5 60 5.0 0.28 | 42 | 0.8 | 33 |0.53 3.6 0.36 
IV 26.0 60 2.3 0.30 | 35 | 0.5 | 28 |0.33 1.9 0.19 


Note: 1) Residual strength parameters shall be reduced by 1/3 for the 
cohesion and 20% for the friction angle; 2) Dilation angle was estimated at 
1/3 of the friction angle; 3) Tensile strength of rock mass is taken as 10% of 
uniaxial compressive strength of rock mass. 
3.2 Shotcrete 

a) Plastic fiber shotcrete: unit weight: 24.0kN/m?; 

b) The Elastic moduli are calculated as formula E.=4700*f.’*0.5 
(ACI 318, section 19.2.2 Modulus of elasticity). 


Table 3.2 Shotcrete strength and elastic modulus 


Compressive strength 


; Elastic modulus E.(MPa) Poisson ratio 
fc. (MPa) 
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25 23500 0.2 
30 25743 0.2 


3.3 Rock Dowels 


Rock dowel type: Cement mortar full length bonded dowel 
1) Yield strength: 500MPa 


2) Compressive/Tensile strength: 545MPa 
3) Elastic modulus: 210,000MPa 
4) Tensile strength of rock bolt steel: 500MPa/1.5-333MPa, 
And the value 1.5 is safety factor. 
Ultimate bearing capacity of rock dowel is calculated from formula 


below: 


F=0.5ndLTur 
Where, 


d =Effective diameter of borehole; 

L = Length of grouted portion of anchor; 

t =Working bond strength; 

Tut= The ultimate bond strength at failure, the ultimate bond stress 
is often taken as 1/10 of the uniaxial compressive strength of the 


rock or grout (whichever is less) (Littlejohn 1977) up to a maximum 
value of 4.2 MPa. 


Table 3.3 The ultimate bearing capacity calculation of rock dowel 


Rock classification Unit II III IV 

UCS MPa 60 60 60 

Rock dowel diameter mm 25 25 25 
Rock dowel area mm? 491 491 491 

Itimat t th of k 42 
Ultimate bond strength of roc kPa 4200 4200 0 

dowel Tut 0 

Effecti i t f k | 
ective diameter of rock dowe um 50 50 50 
borehole d 
Bond strength of unit length MN/ 0.33 0.5mx50x4200 | 0.3 
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/10°= 
F=0.5rdTut m 0.33 i 
Tensile st th of rock | 
ensile strength of rock dowe MPa 333 333 333 
steel 
491x333/109-2 | 0.1 
Tensile peak value of rock dowel MN 0.164 0.164 64 


3.4 Lattice Girder and Reinforcing Bar 
The mechanical parameters of reinforcing bar are shown as below. 


Table 3.4 Mechanical parameter for reinforcing bar 


Yield strength Tension Elastic modulus Poisson 
(MPa) strength(MPa) E.(MPa) ratio 
500 545 210,000 0.3 


3.5 Steel support 


MB150 H-beam type steel rib is used. And the mechanical 
parameters of steel support are shown as below. 


Table 3.5 Mechanical parameter for steel support 


Yield strength Tension Elastic modulus Poisson 
(MPa) strength(MPa) E.(MPa) ratio 
250 400 210,000 0.3 


4 Initial Proposed Tunnel Supports 


4.1 Tunnel Support of Excavation Class 


The initial proposed supports for each rock type are shown as below. 


Table 4.1 Initial proposed supports for each rock type 


Initial proposed supports 
Item Unit 

| II HT IV 
Shotcrete thickness mm Spot 50 80 100 
Shotcrete strength MPa 21 21 21 21 
Rock dowel length m 3.0 3.0 3.0 3.0 
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Rock dowel spacing within each spot spot 2.0 15 


row 
Rock dowel row spacing m spot spot 2.0 1.5 
Steel support/Lattice girder | m spot 


4.2 Safety Factor 
(1) The support safety factor for shotcrete and rock dowel are 1.5. 


(2) The minimum safety factor for wedges analysis is 2.0. 


5 Model Set-up, Boundary and Loading Conditions, Stress 


Initiation 


For simplicity, the model will not extend to the ground surface, and 
the boundary of model is set 10 times tunnel width/height, but to 
apply corresponding stresses at the model’s upper boundary. 

For the case of gravitational loading, only, boundary and loading 
conditions to be set as shown below. This applies for tunnels in the 
UT1-Project that are oriented NNE-SSW (+/- 45°), as NE-SSW is 
generally orientation of the maximum horizontal stress o} in the 
Upper Himalaya, which is in the plane-strain Phase2-model the 
(“irrelevant”) out-of-plane stress. Hence, initial stress conditions in 
the model can be calculated by gravitational loading, or by 
specifying the relative ground surface elevation. 

The finite element model for rock class |, Il, Ill, IV is shown as Figure 
5.1, and the corresponding gravitational loading is shown in the 
Table 5.1. 
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x fixed, 


x fixed, 
y free 


y free 


x fixed, 
y fixed 


Figure 5.1 Finite element model and gravitational loading 


Table 5.1 gravitational loading for surrounding rock 


Overburde 
Tunnel Rock TET : 
Rock class/ btided Bi n Application location of Max. 
Location . gravitation buried depth 
depth | weight 
al loading 
m kN/m? MPa 
| 300 27 8.1 Refer to rock class II location 
ll 300 27 8.1 Whole tunnel 
Junction 300 27 8.1 Class Il 
I 300 26.5 7.95 Whole tunnel 
200 26.5 5.41 Inferred location geologically 
Turni | Hl 
UMS. | sq | "26:5 gg |e 
bay 
200 26 5.20 Refer to rock class III location 
IV 150 26 3.60 Intermediate value 
100 26 2.60 Inferred location geologically 


6 Tunnel Modelling and Results 


6.1 Rock Class I (MAX. Buried Depth 300m) 


For surrounding rock class | 


parameters listed in Section 3, and the calculated displacement and 


plastic zone are shown as below. 


without support, 
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5.10e-004, 


9.10e-004 


Figure 6.1.1 Plastic zone and total displacement without tunnel 
support 
It can be seen that the maximum displacement is 0.5mm in roof and 
0.9mm in invert, and almost no plastic zone. In terms of stability, 


the plastic fiber shotcrete may be required when necessary. 


6.2 Rock Class II (MAX. Buried Depth 300m) 
For rock class Il, according to LDP curve and Characteristic curve 


presented in section 2, it is determined to choose 9m as the max. 
unsupported distance at support installation. The displacement 
without tunnel support, the displacement after support installation, 
plastic zone and axial force of rock dowels, the safety factor of 
plastic fiber shotcrete/ steel rib (including shear and moment) are 


shown as following. 
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“[1.50e-003] 


Figure 6.2.1 Total displacement without tunnel support 


Total 
Displacement 
m 


0.00e+000 
3.00e-004 


6.00e-004 
9.00e-004 
1.20e-003 
1.50e-003 
1.80e-003 
2.10e-003 


2.70e-003 


2.40e-003 
2.70e-003 
3.00e-003 
3.30e-003 


3.60e-003 


Figure 6.2.2 Total displacement after support installation 
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Yielded 
Elements 
[per-cent] 


50 


im 
ES 
EN 
BEN 


Sy 
^ 0.035MN Axial Force - - 
& : 7.987e-004MN Axial Force 


Figure 6.2.3 plastic zone and axial force of rock dowels after support 


installation 


1.5 


Support Element: plastic fiber shotcrete 


Figure 6.2.4 the safety factor of plastic fiber shotcrete 
6.3 Rock Class III (MAX. Buried Depth 300m) 


For rock class Ill, according to LDP curve and Characteristic curve 


presented as following, it is determined to choose 5m as the max. 


unsupported distance at support installation. 
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Characteristic curve of Class III 


Total Displacement ui [m] 
o 
HM 


0 


0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% 


Stress Release Ratio Ai 


Figure 6.3.1 Characteristic curve 


LDP Curve of Class III 
16.00 
14.00 


12.00 


The displacement of vertex (mm) 


-15.0 -10.0 -5.0 ` 0.0 5.0 10.0 15.0 


The distance from tunnel advancing face (m) 


Figure 6.3.2 Longitudinal Displacement Profile curve 
The displacement without tunnel support, the displacement after 
support installation, plastic zone and axial force of rock dowels, the 


safety factor of plastic fiber shotcrete/ steel rib (including shear and 
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moment) are shown as following. 


Total 
Displacement 
m 


0.00e+000 
1.60e-003 
3.20e-003 
4.80e-003 
6.40e-003 


8.00e-003 


9.60e-003 
1.12e-002 1.52e-002 
1.28e-002 
1.44e-002 
1.60e-002 
1.76e-002 


1.92e-002 


Figure 6.3.3 Total displacement without tunnel support 


Total 
Displacement 
m 

0.00e+000 


1.50e-003 
3.00e-003 
4.50e-003 
6.00e-003 
7.50e-003 
9.00e-003 
1.05e-002 
1.20e-002 
1.35e-002 
1.50e-002 
1.65e-002 


1.80e-002 


Figure 6.3.4 Total displacement after support installation 
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Yielded 
Elements 
[per-cent] 
0 


50 


0.116MN Axial Forcel” 


Figure 6.3.5 plastic zone and axial force of rock dowels after support 


installation 


1.5 


v amo i 


Support Element: plastic fiber shotcrete 


Figure 6.3.6 the safety factor of plastic fiber shotcrete 
6.4 Rock Class III (MAX. Buried Depth 200m) 
Some support system is considered: 
(1) The shotcrete thickness is increased to 0.06m; 
Page 28 of 57 


DOOSAN 


UPPER TRISHULI-1 HEP 
(216MW) 


€ 


Cp [s] EB $E 


POWERCHINA 


Surrounding Rock Stability Calculation of Adit No.3 


For rock class Ill, according to LDP curve and Characteristic curve 


presented in section 2, it is determined to choose 5m as the max. 


unsupported distance at support installation. The displacement 


without tunnel support, the displacement after support installation, 


plastic zone and axial force of rock dowels, the safety factor of 


plastic fiber shotcrete/ steel rib (including shear and moment) are 


shown as following. 


Total 


Displacement 


m 


0.00e+ 
1.30e-' 
2.60e- 
3.90e- 
5.20e- 
6.50e- 
7.80e- 
9.10e- 
1.04e- 
1.17e- 
1.30e- 
1.43e- 
1.56e- 


1.18e-002 


1.24e-002 


1.56e-002 


Figure 6.4.1 Total displacement without tunnel support 
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Total 
Displacement 
m 


8.45e-003 


0.00e+000 

1.30e-003 

2.60e-003 

3.902-003 

5.20e-003 

6.50e-003 

7.80e-003 

3.10e-003 d 
1.04e-002 j 
1.17e-002 \ | 


1.30e-002 
1.43e-002 
1.56e-002 


Figure 6.4.2 Total displacement after support installation 


Yielded 
Elements 
[per-cent] 


50 


Figure 6.4.3 plastic zone and axial force of rock dowels after support 


installation 
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1.5 


Support Element: plastic fiber shotcrete 


Figure 6.4.4 the safety factor of plastic fiber shotcrete 


6.5 Rock Class IV (MAX. Buried Depth 100m) 
For rock class IV, According to LDP curve and Characteristic curve 


presented as following, it is determined to choose 3m as the max. 


unsupported distance at support installation. 


Total Displacement ui [m] 


Characteristic curve of Class IV 


1.60E-02 


1.40E-02 


1.20E-02 


1.00E-02 


8.00E-03 
6.00E-03 


4.00E-03 


2.00E-03 


0.00E-- 00 


0% 1096 20% 30% 40% 50% 60% 70% 80% 90% 10096 


Stress Release Ratio Ai 


Figure 6.5.1 Characteristic curve 
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LDP Curve of Class IV 
14.00 


12.00 


10.00 


8.00 


6.00 


4.00 


The displacement of vertex (mm) 


- 0.00 
-15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 


The distance from tunnel advancing face (m) 


Figure 6.5.2 Longitudinal Displacement Profile curve 


The displacement without tunnel support, the displacement after 
support installation, plastic zone and axial force of rock dowels, the 
safety factor of plastic fiber shotcrete/ steel rib (including shear and 


moment) are shown as following. 


Figure 6.5.3 Total displacement without tunnel support 


Page 32 of 57 


DOOSAN UPPER TRISHULI-1 HEP *$ CO EER 2 


(216MW) mat Bay 


Surrounding Rock Stability Calculation of Adit No.3 


"7.80e-003] 


Figure 6.5.4 Total displacement after support installation 


Yielded 

Elements 

[per-cent] 
0 


-0.005MN Axial Force 


50 


100 


Figure 6.5.5 plastic zone and axial force of rock dowels after support 
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installation 


12 
1.5 


Fhe 


Qe 


06 05 04 03 02 01 00 Of 02 03 04 05 06 OF 


Support Element: plastic fiber shotcrete 


Figure 6.5.6 the safety factor of plastic fiber shotcrete 


6.6 Rock Class IV (MAX. Buried Depth 150m) 
For rock class IV, According to LDP curve and Characteristic curve 
presented as following, it is determined to choose 3m as the max. 


unsupported distance at support installation. 


Characteristic curve of Class IV 
4.00E-02 
3.50E-02 
3.00E-02 
2.50E-02 
2.00E-02 
1.50E-02 
1.00E-02 
5.00E-03 
0.00E--00 


Total Displacement ui [m] 


0% 10% 20% 30% 4096 50% 60% 70% 80% 90% 10096 


Stress Release Ratio Ai 
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Figure 6.6.1 Characteristic curve 


LDP Curve of Class IV 
35.00 


30.00 


25.00 


20.00 


15.00 


The displacement of vertex (mm) 


-15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 


The distance from tunnel advancing face (m) 


Figure 6.6.2 Longitudinal Displacement Profile curve 
Some support system is considered: 
(1) The shotcrete thickness is increased to 0.12m; 
(2) The shotcrete strength is increased to 30MPa; 
(3) The spacing of Steel support is taken as 1m; 
(4) The spacing of rock dowel is decreased to 1m; 
The displacement without tunnel support, the displacement after 
support installation, plastic zone and axial force of rock dowels, the 
safety factor of plastic fiber shotcrete/ steel rib (including shear and 


moment) are shown as following. 
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Total 3.37e-002 
Displacement 
m 


0.00e+000 
4.00e-003 
8.00e-003 
1.20e-002 
1.60e-002 
2.00e-002 
2.40e-002 


2.80e-002 7 3.00e-002 
3.20e-002 


3.60e-002 
4.00e-002 
4.40e-002 


4.80e-002 


Figure 6.6.3 Total displacement without tunnel support 


Total 
Displacement 
m 

0.00e+000 


3.00e-003 
6.00e-003 
9.00e-003 
1.20e-002 
1.50e-002 
1.80e-002 


2.10e-002 
2.40e-002 
2.70e-002 
3.00e-002 
3.30e-002 


3.60e-002 


Figure 6.6.4 Total displacement after support installation 


Page 36 of 57 


DOOSAN UPPER TRISHULI-1 HEP 9% CEB c3 


(216MW) mat Bai 


Surrounding Rock Stability Calculation of Adit No.3 


Yielded 

Elements 

[per-cent] 
0 


50 


Figure 6.6.5 plastic zone and axial force of rock dowels after support 


installation 


07 06 05 04 03 02 01 00 01 02 03 04 05 06 07 08 4 
Shear force Q - MN Moment M - MNm 


Support Element: plastic fiber shotcrete 


Figure 6.6.6 the safety factor of plastic fiber shotcrete 
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6.7 Rock Class IV (MAX. Buried Depth 200m) 

For rock class IV, according to LDP curve and Characteristic curve 
presented in section 2, it is determined to choose 3m as the max. 
unsupported distance at support installation. 

Some support system is considered: 

(1) The shotcrete thickness is increased to 0.16m; 

(2) The shotcrete strength is increased to 30MPa; 

(3) The spacing of Steel support is taken as 1m; 

(4) The spacing of rock dowel is decreased to 1m; 

(5) The length of rock dowel is increased to 3.5m. 

The displacement without tunnel support, the displacement after 
support installation, plastic zone and axial force of rock dowels, the 
safety factor of plastic fiber shotcrete/ steel rib (including shear and 


moment) are shown as following. 


(6.41e-002 


5.95e-002 


7.70e-002 


Figure 6.7.1 Total displacement without tunnel support 
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Total 
Displacement 
m 

0.00e+000 


6.00e-003 
1.20e-002 
1.80e-002 
2.40e-002 
3.00e-002 
3.60e-002 
4.20e-002 
4.80e-002 
5.40e-002 
6.00e-002 
6.60e-002 
7.20e-002 


3.60e-002 


6.30e-002 


Figure 6.7.2 Total displacement after support installation 
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Figure 6.7.3 plastic zone and axial force of rock dowels after support 


installation 
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Legend 
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Support Element: plastic fiber shotcrete 


Figure 6.7.4 the safety factor of plastic fiber shotcrete 
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beam): W150x18 Teagan 


EN 


Support Element: steel rib 


Figure 6.7.5 the safety factor of steel rib 
6.8 Turning Bay (MAX. Buried Depth 150m) 
For turning bay, it is placed in rock class III (Its position is artificially 
controlled), according to LDP curve, the Characteristic curve 
presented as following, it is determined to choose 5m as the max. 


unsupported distance at support installation. 


Characteristic curve of Class III 


Total Displacement ui [m] 


0d 
0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 10096 


Stress Release Ratio Ai 


Figure 6.8.1 Characteristic curve 
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LDP Curve of Class Ill 
3-50 


3.00 


2.50 


The displacement of vertex (mm) 


0.00 
-15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 


The distance from tunnel advancing face (m) 


Figure 6.8.2 Longitudinal Displacement Profile curve 
The displacement without tunnel support, the displacement after 
support installation, plastic zone and axial force of rock dowels, the 
safety factor of plastic fiber shotcrete/ steel rib (including shear and 


moment) are shown as following. 


Total 
Displacement 
m 

0.00e+000 
4.00e-004 3.28e-003 
8.00e-004 
1.20e-003 
1.60e-003 
2.00e-003 
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2.80e-003 
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3.20e-003 
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4.40e-003 


4.80e-003 4 40e-003 
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Figure 6.8.3 Total displacement without tunnel support 
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Displacement 
m 
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8.00e-004 
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> paden 2.80e-003] 2.80e-003 
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Figure 6.8.4 Total displacement after support installation 
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Figure 6.8.5 plastic zone and axial force of rock dowels after support 


installation 
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Legend 


1.2 
1.5 


Support Element: plastic fiber shotcrete 


Figure 6.8.6 the safety factor of plastic fiber shotcrete 


6.9 Junction (MAX. Buried Depth 300m) 

The junction is inferred as rock class Il, according to LDP curve, the 
Characteristic curve presented as following, it is determined to 
choose 9m as the max. unsupported distance at support 


installation. 


Characteristic curve of Class Il 


Total Displacement ui [m] 
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Figure 6.9.1 Characteristic curve 
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LDP Curve of Class II 
3-00 


2.50 


2.00 


1.50 


1.00 


The displacement of vertex (mm) 


0.00 
-15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 


The distance from tunnel advancing face (m) 


Figure 6.9.2 Longitudinal Displacement Profile curve 
The displacement without tunnel support, the displacement after 
support installation, plastic zone and axial force of rock dowels, the 
safety factor of plastic fiber shotcrete/ steel rib (including shear and 


moment) are shown as following. 
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Figure 6.9.3 Total displacement without tunnel support 
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Figure 6.9.4 Total displacement after support installation 


Yielded 
Elements 
[per-cent] 


50 


Figure 6.9.5 plastic zone and axial force of rock dowels after support 


installation 
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Support Element: plastic fiber shotcrete 


Figure 6.9.6 the safety factor of plastic fiber shotcrete 


7 Empirical support categories 


7.1 Excavation support ratio (ESR) 
According to the handbook of the Q-system, a low ESR value 
indicates the need for a high level of safety while higher ESR values 
indicate that a lower level of safety will be acceptable. 
Requirements and building traditions in each country may lead to 
other ESR-values than those given in table 7.1. 

Table 7.1 ESR-values 


Page 47 of 57 


DOOSAN UPPER TRISHULI-1 HEP 4% COE zi 


(216MW) mat Bay 


Surrounding Rock Stability Calculation of Adit No.3 


7 Type of excavation ESR 


Temporary mine openings, etc. 


Vertical shafts*: i) circular sections 
ii) rectangular/square section 


* Dependant of purpose. May be lower than given values. 


C Permanent mine openings, water tunnels for hydro power (exclude high pressure 16 
penstocks), water supply tunnels, pilot tunnels, drifts and headings for large openings. ! 


- Minor road and railway tunnels, surge chambers, access tunnels, sewage tunnels, etc. ez 


Power houses, storage rooms, water treatment plants, major road and railway tunnels, 
civil defence chambers, portals, intersections, etc. 


Underground nuclear power stations, railways stations, sports and public facilitates, 
factories, etc. 


Very important caverns and underground openings with a long lifetime, « 100 years, or 
without access for maintenance. 


As for construction adit tunnel, the ESR=1.3. 

7.2 Excavation support categories 

According to rock parameters for tunnels, rock class | Q>40, rock 
class Il Q>10 & Q<40, rock class Ill Q>4 & Q<10, rock class IV Q>1 
& Q<4. 

The span of tunnel is about 5m. 

The equivalent dimension=Span or height in meter/ESR=5/1.354m. 


As for class I~IV, the empirical supports are shown as follows. 
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Figure 7.2 Empirical support for rock class II 
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Figure 7.4 Empirical support for rock class IV 
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Figure 7.5 Empirical support for junction 
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Figure 7.6 Empirical support for turning bay 
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Table 7.2 Empirical proposed supports 


Empirical proposed supports 


Item Unit : : 
| II unction Hl Turning IV 
Shotcrete thickness mm | spot | spot 50 50 50 60 
Rock dowel length m | spot | spot | 2.6 2.4 2.7 2.4 
Rack dawel TOW m | spot | spot 2.0-3.0 |-97?: | 1.6~2.0 | 1.6 
spacing 0 


8 Wedge Analysis 


8.1 Discontinuity Input Data 

According to the site investigation, there are three group joints may 
distribute in the adit tunnel, as is shown in Table 8.1. The 
stereographic projection of these joints in rock mass in shown in 
Figure 8.1~8.4. The Mechanical parameter of joints is shown in Table 
8.2. 


Table 8.1 Occurrence of joints in adit tunnel 


Joint NO. Dip Dip Direction Tunnel Direction 
Jl (foliation) 25° 333° 
J2 65° 160° 
N29°W 
J3 75° 120° 
J4 60° 65° 


Table 8.2 Mechanical parameter of joints in the adit tunnel 


JCS C ® 
JONIT TYPE JRC 
[L]MPa[] [L]MPa[] EE 
Foliation(J1) 6 70 0.20 35.0 
lithoclasts 
cemented 8 50 0.15 28.8 
(020304) 
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Figure 8.1 Stereographic projection of the joins (J1/J2/J3) in the adit 


tunnel 
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Figure 8.2 Stereographic projection of the joins (J1/J2/J4) in the adit 


tunnel 


Page 54 of 57 


UPPER TRISHULI-1 HEP 
minimum (216MW) <> S228 


Surrounding Rock Stability Calculation of Adit No.3 


Figure 8.3 Stereographic projection of the joins (J1/J3/J4) in the adit 


tunnel 
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Figure 8.4 Stereographic projection of the joins (J2/J3/J4) in the adit 


tunnel 


8.2 Wedge Stability Analysis 
According to the wedge stability analysis, potential wedges may 
occur in the excavation periods, which are shown in Table 8.3. 


Table 8.3 Instability Block Search Statistics of adit tunnel 


Fissur | Block | Maximu | Safe 
Developm | Instabili Graphical 
S. | eset | volu | m block | ty 
. ent ty presentation of 
N. | numb me height | fact " 
: position mode block 
er (m?) (m) or 
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J1J2/ 0.82 | Lowerleft | |... SS 
1 291 1. | i 
J3 d o 4. | sidewall je \ | 
Upper y 
2 JA | 43901] 2.26 | 973 right sliding | 
J3 2 | 
side wall C 
PF» 
3 TUE 0.154 1.16 0:09 Crown fall | | 
J3 8 | | 
Lower y 
J1J2/ | 12.81 320| . May be 
4 8.92 ht sid 
J4 4 6 HS SIS stable 
wall 
J1/]2/ 0.36 | Upper left -— 
.11 12.2 | 
> [j4 ii 9 | a l sidewall |e 
After supporting, shotcrete | 4.15 May be ~ 
50mm 6 stable 
A 
L be 
6 re 0.087 | 2.51 0 Crown fall 
JJ3/ | 37.20 Power 
7 j4 b 10.45 | 89.5 | rightside | sliding 
wall 
After supporting, shotcrete May be Pp. 
190 ^ ^. 
50mm stable 
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Juja/ o.74 | “PPE 7 p. 

8 0.1 0.47 right sliding FA 
J4 9 . 
side wall <1 
Jl/J3/ | 69.55 Upper left | |... ~ 
9 j4 ^ 13.01 0.7 Sides wali sliding N 
After supporting, shotcrete May be VG > 
2.99 
50mm stable 
1 
10 "d 0.015 1.98 0 Crown fall 
J2/)3/ Lower left m à 
11 .224 1.71 | | 
JA : "adepan | 
Upper 
2 i 
12 uo 0.306 1.99 yo right sliding 
side wall 
13 d 6.895 2.88 ids Crown fall 
J4 6 
After supporting, shotcrete 4.26 May be k 
50mm 5 stable X | 


It could be seen that, in all the possible joints combination except 


No.5/7/9/13, the volume of the unstable block is about 0.1—1.3m^, 


and the height of the block is about 0.4m~2.5m. In the tunnel, as 


for class | & 2, spot shotcrete can restrain the trend of sliding and 
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fall. As for other classes, the pattern bolting can control the sliding 
and fall. During the excavation process, a reasonable layered 
excavation sequence shall be used and a reasonable excavation 
height to avoid large potential unstable blocks. 

As for the NO.5/7/9/13 joint combination, after spot and pattern 
supports, the minimum safety factor of potential wedges without 
supports in adit tunnel is bigger than the minimum safety factor of 
2.0. 


9 Conclusions and Recommendations 


9.1 Conclusions 

As per the calculation fore adit tunnel above, it can be concluded as 
follows: 

(1) For rock class I, It can be seen that the maximum displacement 
is less than 0.1% of tunnel height (0.001x5600=5.6mm), the 
systematic support is not required. In the actual excavation process, 
spot shotcrete can be carried out according to the exposed rock. 

(2) For rock class Il, the tunnel roof occur plastic zone. In case of fall- 
block, the spot rock dowels many be required in the tunnel roof. And 
the thickness of shotcrete should be 50mm thick. 

(3) For rock class IV, while the maximum buried depth is no greater 
than 100m, the initial proposed supports meets the requirement of 
tunnel stability; while the maximum buried depth is greater than 
100m and less than 150m, the thickness of shotcrete should be 
improved to 120mm and the strength of shotcrete should be 
improved to 30MPa and the rock dowel spacing is decrease to 1m; 


while the maximum buried depth is greater than 150m and less 
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than 200m, more measures should be provided, e.g.: The shotcrete 
thickness is increased to 160mm, the shotcrete strength is 
increased to 30MPa, the spacing of Steel support and rock dowel is 
decreased to 1m. 

(4) For wedge analysis, after shotcrete, the minimum safety factor 
of potential wedges without supports in adit tunnel is bigger than 
the minimum safety factor. However, the drainage system in 
construction periods is extremely important in the stability 
calculation of wedges, it should be paid more attention in the 


excavation process. 


9.2 Recommendations 
According to the calculation above, the recommended supports for 
each rock type are summarized in the Table 9-1~9.2. 


Table 9.1 Recommended supports for each rock type 


Rock class | Max. Shotcret | Rock Rock Steel 
buried e dowel dowel support/Lattice 
depth(m) thicknes | length row girder spacing 
S (m) spacing( | (m) 
(mm) m) 
I| SC 300 spot 3.0 spot / 
SE NA / / / / 
SP 300 spot | 3.0 spot |/ 
Il | SC 300 spot 3.0 spot / 
SE NA / / / / 
SP 300 50 3.0 2.5(ro |/ 
of) 
Il | SCL 300 80 3.0 2.0 / 
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| | SC2 200 60 3.0 2.0 / 
0 

SP 70~300 | 70 3.0 2.0 Z 
| | SC1 200 160 3.5 1.0 1.0 
V | SC2 150 120 3.0 1.0 / 

SC3 100 100 3.0 1.5 / 

SE NA 60 2.4 1.6 / 

SP 70—200 | 120 3.0 1.5 spot 


Note: SE denotes empirical support, SC denotes calculating support, 
SP denotes proposed support. 


Table 9.2 Recommended supports for junction & turning 


Junction(ll) Turning bay(lll) 
Item Unit 
SC SE SP |SC SE SP 
Effective equivalent span m 7.2 8.0 

Shotcrete thickness mm|80| 50 [80 | 10| 50 | 22 

n n 
Rock dowel length m 13.0 2.6 3.0 | 3 2.7 3.0 
Rock dowel row spacing m |2.0 x 2.0 |2.0 ui 2.0 
Steel support/Lattice girder spo Spo 

l m | / / / / 

spacing t t 

Stability Yes | Yes |Yes Yes | Yes Yes 
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